Cores in Dwarf Galaxies from Dark Matter with a Yukawa Potential 
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We show that cold dark matter particles interacting through a Yukawa potential could naturally 
explain the recently observed cores in dwarf galaxies without affecting the dynamics of objects with 
a much larger velocity dispersion, such as clusters of galaxies. The velocity dependence of the 
associated cross-section as well as the possible exothermic nature of the interaction alleviates earlier 
concerns about strongly interacting dark matter. Dark matter evaporation in low-mass objects 
might explain the observed deficit of satellite galaxies in the Milky Way halo and have important 
implications for the first galaxies and reionization. 



PACS numbers: 95.35+d,98.80-k,98.62-g 

Introduction. The collisionless cold dark matter 
(CDM) model has been highly successful in account- 
ing for the gravitational growth of density perturba- 
tions from their small observed amplitude at early cosmic 
times (as imprinted on the cosmic microwave background 
anisotropies [I ) to the present-day structure of the Uni- 
verse on large scales. However, it is far from clear that 
the predictions of this model arc valid on small scales. 

New data on low mass galaxies indicate that their 
dark matter distribution has a core [2J, in contrast to 
the cusped profile expected from collisionless CDM sim- 
ulations [3]. The mean value of the inner logarithmic 
slope of the mass density profile in seven dwarf galaxies 
within the THINGS survey is observed to be -0.29±0.07 
[3], much shallower than the expected slope of ~ — 1 
from pure CDM simulations. Moreover, the dynamics of 
dwarf spheroidal galaxies, such as Fornax [5], Ursa-Minor 
[Hj, and Sculptor [7], whose luminosities and dynamical 
masses are smaller by 2-3 orders of magnitude than the 
THINGS galaxies, indicates a characteristic core density 
of - O.1±O.O5M pc~ 3 = (7 ±4) x 1CT 24 g cm- 3 . Since 
these dwarf spheroidals are dominated by dark matter 
throughout, it is challenging to explain their inferred 
cores by the gravitational interaction of the dark mat- 
ter with the baryons [8]. Although it is conceivable that 
powerful gas outflows from an early baryon-dominated 
nucleus would reduce the central dark matter density 
in luminous galaxies [HI [TD] , the formation of a massive 
baryonic nucleus would initially compress the CDM [TT] 
and exacerbate the discrepancy that needs to be resolved 
[5], and also potentially violate the observed low lumi- 
nosities from dwarf galaxies at higher redshifts [T2J [T3] . 
High-rcdshift observations of dwarf galaxies must find ev- 
idence for the required strong feedback phase, or else an 
alternative process is at play. Some recent simulations 
that include feedback do not observe the appearance of 
cores within the lowest luminosity galaxies [14] . 

To alleviate early signs of the above discrep- 
ancy, Spergel & Steinhardt [TS] adopted the Strongly- 



Interacting Dark Matter (SIDM) model pMIH] in which 
the dark matter has a large cross-section for self inter- 
action. It was expected that if dark matter scatters in 
the cores of galaxies, then it might resemble a fluid with 
a flatter central density profile. The SIDM proposal fell 
out of favor because: (i) gravitational lensing and X-ray 
data indicate that the cores of clusters of galaxies are 
dense and ellipsoidal, whereas SIDM would predict them 
to be shallow and spherical [HO [5D]; (ii) relaxation of 
the core ultimately generates an even denser nucleus af- 
ter many collision times due to the gravothermal catas- 
trophe, familiar from core collapse of globular clusters 
|21) . although this evolution might take more than the 
Hubble time; (Hi) dwarf galaxies would be expected to 
evaporate when interacting with the higher velocity par- 
ticles of their host halo; and (iv) theoretical biases sug- 
gested that the required cross section was incompatible 
with popular models of Weakly-Interacting Massive Par- 
ticles (WIMPs). 

Recently, there has been growing interest in the pos- 
sibility that WIMPs exhibit "dark forces" as a means 
to address a wide range of anomalies [22] . In particu- 
lar, it was realized that a new force carrier <j> (scalar or 
vector) might naturally mediate a long-range interaction 
on the scale of the de-Broglie wavelength of the WIMPs, 
leading to a self-interaction cross-section for scattering 
that is much greater than for WIMP annihilation. The 
studied forces have a variety of scales in them, from the 
screening scale set by the mass of the carrier particle 
to the non-perturbative scale set by its coupling, adm,p- 
Moreover, these forces are naturally accompanied by new 
energy states. Up- and down-scattering processes (which 
are endo- and exo-thermic, respectively) naturally intro- 
duce yet another scale into the problem, set by the en- 
ergy splitting S between the states. While similar models 
have been previously considered [22] , our regime of inter- 
est was assumed to be unsuitable due to WIMP capture 
and WIMP annihilation. As it turns out, these processes 
are generally weaker than scattering and so this range of 
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parameters remains open. Nonetheless, the absence of 
dramatic departures from CDM predictions has allowed 
important constraints to be placed [24] [25] . 

In this Letter, we examine the possible existence of a 
dark force from a different perspective. Rather than limit 
its allowed range of parameters based on observations, 
we show that it can ameliorate tensions in astrophysi- 
cal data. In particular, we find that a Yukawa force in 
dark matter scattering would naturally produce cores in 
dwarf galaxies while avoiding the myriad constraints on 
SIDM which arise in systems with a much larger veloc- 
ity dispersion, such as clusters of galaxies. The specific 
velocity dependence of the interaction cross-section, as 
well as the possible exothermic nature of the interaction, 
alleviate earlier concerns about the SIDM model. To dis- 
tinguish from previous approaches with a constant cross 
section or a simple power law velocity dependence, we 
label this scenario as Yukawa-Potential Interacting Dark 
Matter (YIDM). 

Dark Forces. The mediator of the force <f> could be 
either a scalar or a vector, as magnetic-type interactions 
are negligible. The force could couple to standard model 
fields through kinetic mixing with the photon, or through 
mass mixing with the Higgs boson. Constraints on the 
presence of such a force come from a wide range of pro- 
cesses [27], but ample parameter space remains for 
a small mixing angle, e ~ 10~ 3 . New searches are un- 
derway to find precisely such a force carrier at ~ GeV 
energy experiments [28] . 

Scattering through a massive mediator is equivalent to 
having a Yukawa potential. The elastic scattering prob- 
lem is then analogous to the screened Coulomb scatter- 
ing in a plasma which is well fit by a cross-section 
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where j3 = irv^/v 2 — 2ad'm < j > / (m x v 2 ) , and v is the rela- 
tive velocity of the particles. We use angular brackets to 
denote that this is the momentum-transfer weighted cross 
section. Here, v a is the velocity at which the momentum- 
weighted scattering rate (o~v) peaks at a cross section 
value of (T m ax = 22.7/m^. The above expression can be 
approximately generalized to the inelastic case by sub- 
stituting — > yjm x 8 for the characteristic minimum 
momentum transfer when < y/m x S (see discussion 
in [30 J. This expression is derived using classical physics, 
and thus, it is important to note what quantum effects 
can come into play. In cases where the de Broglie wave- 
length is longer than the Compton wavelength of the 
force 777^ , the quantum calculation should be consid- 
ered for quantitative results. Nonetheless, the same qual- 
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FIG. 1: Dependence of the self- interaction cross-section (a) on 
the relative velocity (v) for dark matter interacting through a 
Yukawa potential. The normalizations of a and v are set by 
free parameters in the underlying Lagrangian (see Appendix) , 
and we show two possible curves peaking at v a = 1 km s~ 
and = fOO km s _1 [blue, solid and purple, dashed, respec- 
tively) . 

itative features should remain: the cross section should 
saturate at low velocities near a ~ "i^ 2 , and at high 
velocities, where the classical approximation is valid, it 
should fall rapidly. 

Figure 1 depicts the velocity dependence of the elas- 
tic cross-section in Eq. ([I]). Interestingly, the scattering 
rate is nearly constant at low velocities, peaks at a ve- 
locity v a , and declines sharply at v > v a , allowing it to 
introduce cores in dwarf galaxies where the velocity dis- 
persion is low (v ~ 10 km s _1 ) but not in clusters of 
galaxies where the characteristic velocities are larger by 
two orders of magnitude (v ~ I0 3 km s _1 ). The nor- 
malizations of the cross-section and velocity are deter- 
mined by free parameters in the interaction Lagrangian 
(see Appendix) , with the Compton wavelength of the in- 
teraction setting the relevant spatial scale. We show two 
possible values of the peak velocity, one that would pro- 
duce cores only in dwarf galaxies {v a = 10 km s _1 ), and 
another that would produce cores in more massive galax- 
ies (v a = I0 2 km s _1 ) as implied by data on low surface 
brightness galaxies 31 . At any given halo mass, we ex- 
pect scatter in the core properties of individual halos, 
due to variations in their age and assembly history. 

Having one collision per Hubble time at the character- 
istic core density of dwarf galaxies ~ 0.1M Q pc -3 , trans- 
lates to the condition (777 x /IOGeV)(777 /IOOMeV) 2 ~ 1 
(see Appendix). An order of magnitude larger cross- 
sections are also allowed by the data. Figure 2 shows 
the allowed parameter ranges [35] that would naturally 
explain the dark matter distribution in observed astro- 
physical objects. We find that even though collisions 
shape the central profiles of dwarf galaxies, the standard 
collisionless treatment still provides an excellent approx- 
imation for the dark matter dynamics in X-ray clusters. 
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FIG. 2: Astrophysical constraints on the normalization of the 
self-interaction cross-section (a max ) as a function of the ve- 
locity at which the peak collision rate is obtained (« CT ) in Fig. 
1. The red solid line is normalized to have (a) /m v « 

\ / max i a. 

6 x 10~ cm /GeV at tMwarf ~ 10 km s~ , which should be 
regarded as the minimum interaction necessary to flatten the 
cores of dwarf galaxies. Additional lines indicate upper lim- 
its on the cross-section based on astrophysical considerations: 
X-ray cluster ellipticity (blue, dashed), limiting (<J max /m x ) ~ 
4 x 1(T 26 cm 2 /GeV at v ~ 10 3 km s" 1 ; destruction of dwarf 
sub-halos through collisions with high velocity particles from 
a larger parent halo in which these dwarfs are embedded 

(green, dotted), limiting (<r max /ni x ) ~ 5 x 10~ 25 cm 2 /GeV 
|32| at v ~ 200km s~ ; and requiring the number of scatters 
in dwarfs to be less than ~ 10 2 during the age of the Universe 
to avoid the gravothermal catastrophe (purple, dash-dotted) . 
Related limits are summarized in [25] . 



Exothermic interactions. The presence of excited 
states related to a new force has important implications 
for the properties of dark matter. In particular, upscat- 
tering ("inelastic dark matter" (iDM) [33]) or downscat- 
tering [3"4"H58] off nuclei can dramatically affect direct de- 
tection experiments. Dark matter self scattering into an 
excited state ("eXciting Dark Matter" (XDM) [39]) has 
been invoked to explain the excess 511 keV flux observed 
by INTEGRAL [31 [3JHH1 S3 S2] • 

We focus here on the response of dwarf galaxies to 
the presence of excited states, YIDM*, which we assume 
are copiously present from the early universe. In par- 
ticular, the release of kinetic energy in YIDM* collisions 
would help to evade the gravothermal catastrophe [3T] , in 
just the same way that the energy released by primordial 
stellar binaries weakens core collapse in globular clusters 
[43]. 

If excited states exist, then a major fraction of the 
CDM might be excited when the WIMPS decouple ther- 
mally in the early Universe. This excitation could be 



stable on cosmological times in models where the dark 
force mixes with electromagnet ism [35 135] . The scatter- 
ing process can decouple in the early universe at temper- 
atures above the splitting, leaving essentially equal abun- 
dances of all states of the dark matter. Alternatively, the 
excited states could be produced non-thermally. The for- 
mer scenario tends to require light (m x ~ MeV) particles, 
while the latter would be more naturally weak-scale. Af- 
ter decoupling, the early CDM dynamics is identical to 
the standard collisionless model until dwarf galaxies form 
and the crossing of dark matter streams at low v occurs in 
their cores, giving rise to self-interactions on a timescale 
shorter than the age of the Universe. 

We denote the velocity imparted to an initially slowly- 
moving WIMPs (v < Vcrit) upon scattering by u cr ;t = 
y S/m. Since the gravitational potential of a dwarf 
galaxy halo is shallow, sufficiently exothermic collisions 
could eject colliding particles out of the halo if w cr it ex- 
ceeds the local escape velocity. If XX* ~> XX scatterings 
can occur (see Appendix), then the halo core will lose 
particles until it reaches a density such that the interac- 
tion time is comparable to the age of the halo. Requiring 
that the final core particles will interact only once over 
the current age of the Universe yields a final mass density 
of dark matter, 



Px 



0.1' 



M Q / 0-/6x10" 



pc 



m x /GcV 



10 km s 



•(2) 



The profile of p x will then be set by the velocity depen- 
dence of (av) and the gravitational potential in a steady 
state. Below we show that for an exothermic interaction 
(av) is constant, leading naturally to a constant density 
core in dwarf galaxies. This model predicts that dwarf 
galaxies of a similar age should have a similar core den- 
sity in agreement with interpretations of data for the 
nearest dwarf galaxies [44] . 

If the characteristic scattering velocity is much higher 
than i> cr it, the process is essentially elastic and one can 
employ Eq. ([!]). On the other hand, if the process is 
highly exothermic then the scattering rate is given sim- 
ply by av = 27ri> cr i t a^/m^ for m x S < (assuming the 
scattering process is still perturbative, i.e. f3 < 1). The 
resulting velocity-independent av would naturally pro- 
duce cores with a flat density profile in dwarf galaxies. 

The density flattening in dwarf galaxies does not im- 
ply the same upper limit on the dark matter density in 
all its cusps. For massive halos, the release of excess 
kinetic energy by collisions has a marginal significance, 
since it only perturbs the low-velocity tail of the CDM 
distribution function and adds a negligible energy at high 
relative velocities where the majority of particles have a 
low interaction rate anyway. 

The evaporation of exothermic YIDM* from dwarf 
galaxies with a gravitational binding energy below the 
energy released in collisions could potentially account for 
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the deficit in the observed abundance of dwarf galaxies 
relative to theoretical CDM expectations [45]. Numer- 
ical simulations are necessary to reliably quantify this 
important effect. Also, dark matter halos which accrete 
cold gas at early times but evaporate at late time might 
leave behind a star cluster with little dark matter. If so, 
some old globular clusters [IS] might be the sought-after 
remnants of the missing dwarf galaxies in the Milky Way 
halo. 

Evolution with redshift. The primordial density per- 
turbations are modified by WIMP scatterings, but for 
m x ~ GeV this modification ends well before observ- 
able modes enter the horizon, leaving the standard nearly 
scale-invariant power-spectrum of density perturbations. 
The imprint of collisions on the density profile of halos is 
expected to evolve with redshift, because at earlier cos- 
mic times halos are denser and younger. A halo of mass 
M collapsing at a redshift z > 1 has a characteristic 
virial radius [IT], 

a corresponding circular velocity, 

and an age limit of ~ 0.5 Gyr[(l + z)/10]~ 3 / 2 . It would 
be particularly interesting to explore the formation of the 
first galaxies using numerical simulations of the YIDM 
model. Deviations from the standard CDM predictions 
could be tested by upcoming galaxy surveys or 21-cm ob- 
servations of the high-redshift Universe j47j [48] . For ex- 
ample, the scaling a oc v~ A might not be allowed to con- 
tinue to arbitrarily low velocities, as this would have de- 
layed reionization beyond observational constraints [49) . 

Discussion. We have shown that a Yukawa potential 
interaction of dark matter can explain the recent data on 
dark matter cores in dwarf galaxies [H [7J , while evading 
the many constraints previously considered for SIDM. 
The new ingredients of the YIDM model involve the non- 
trivial velocity dependence of the scattering rate and the 
possibility of exothermic interactions. 

The velocity dependence of the cross-section does not 
have a simple power-law form, as invoked previously 
[2"01 . The presence of a plateau with a sharp cut-off 
in the scattering rate of YIDM allows the interaction to 
be effective for dwarf galaxies while being entirely sup- 
pressed at high velocities relevant for cluster cores and 
the evaporation of small sub-halos within bigger halos. 
At the same time, it does not rise indefinitely at low ve- 
locities and thus avoid other concerns |32j . 

Excited states naturally accompany a dark force and 
introduce a qualitatively new ingredient: the possibility 
of energy release. In deep gravitational potentials the 



scatterings will be similar to the standard elastic case, 
whereas in shallow gravitational potentials a single scat- 
ter can eject the colliding dark matter particles from the 
halo. Exothermic interactions drive dwarf halos to con- 
stant density cores. Moreover, for the earliest forming 
halos with the highest dark matter density, it is possible 
that the dark matter would evaporate at later times, leav- 
ing behind baryonically-dominated stellar clusters such 
as the oldest globular clusters. 

Unlike SIDM with a single parameter of a/m x , YIDM 
has a broader parameter space. Out of the possible com- 
binations of the underlying particle physics parameters 
(a,5 } m x ,m^), only three affect the properties of halos: 
the cross-section cr max and velocity v a at which the scat- 
tering rate (av) peaks; and in the YIDM* scenario there 
is the additional parameter of the velocity imparted to a 
particle undergoing an exothermic scattering, w cr it- 

Since halos form hierarchically and the coarse-grained 
phase-space density can only decrease through their 
mergers (Liouville's theorem), the development of a core 
in dwarf galaxies would trim the central cusp in big- 
ger halos as well. Numerical simulations are required 
to quantify the average behavior as well as the scatter 
in the core properties of halos as a function of mass and 
redshift. 

Finally, we note that the similarity between the re- 
quired self-interaction cross-section per unit mass of the 
dark matter and baryons may indicate a deep underlying 
relationship between these components. 

Appendix: Particle Physics Implications 

Although a thorough investigation of the particle 
physics parameter space is beyond our scope here, it is 
worthwhile to at least briefly consider the range of param- 
eter space available to us. As mentioned, the mediator 
of the force 4> could be either a scalar or a vector, since 
magnetic-type interactions are negligible. Simple pertur- 
bative realizations of dark forces involve models with an 
additional boson 0^ of a small mass ~ GeV. We 
can assume that the WIMPs are charged under a new 
C^(l)dark symmetry and allow the possibility of a small 
splitting between their Majorana states, leading to the 
Lagrangian, 

C = xPx + ^F^F^ + eF^F^+ml^r 

+ ™ X XX + S XX- (5) 

This can be trivially generalized to a non-Abelian model 
with multiple excited states [HI [50] , which induces split- 
tings between the states radiatively at order o^m^, where 
ad is the fine structure constant of the dark force. Simi- 
larly, the force can arise through a scalar exchange with 
the additional terms 

C D yxx^ + P W + m$0> (6) 
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Here, the effective coupling constant ay — y 2 /Air. Elastic 
scatterings can arise for both scalar and vector, while 
inelastic scatterings are most natural in the gauged case. 
We begin by focusing on the elastic scattering scenario. 

The natural mass scale that determines the maximum 
momentum transfer cross section is m§. For large val- 
ues of 0.1 ~ /3 ~ 10 3 (above which the cross section in 
galaxy clusters is ~ 10~ 3 that in dwarf galaxies), Eq. 
Q yields (a) w (167r/m|)/? 2 /(2 + 3/3 165 ). Requiring 
a cross section/mass of 6 x 10~ 25 cm 2 GeV -1 then im- 
plies m x m 2 = 3 x l(T 2 ^ 2 /(2 + 3/3 165 ) GcV 3 . Thus, for 
m x = 10 (10 3 ) GeV, one has m < 10 2 (10) MeV. For 
standard WIMPs interacting through weak scale forces, 
such phenomena do not occur; however, for dark matter 
augmented by light dark forces it is possible. For lighter 
(~ GeV) WIMPs, a quantum calculation should be used 
to generate the precise quantitative results. The relevant 
mass range for <f> is not surprising in light of the relatively 
weak constraints from galactic dynamics found by recent 
studies [25 | I25 ] . 

Adopting /3 = /3dwarf in a typical dwarf 
galaxy implies through the above relation, 
a d « 10- 7 (m x /GeV) 3 / 2 ^3/3^/2 + 1. The de- 
pendence on /?dwarf and m x allows a wide range of 
values, 10 -5 ~ ay ~ 1, with the low end for /3dwarf ~ 1 
and m x ~ lOGeV, and the high end for /?dwarf ~ 10 3 
and m x ~ 10 3 GeV. 

Although ad can in principle have any value in the 
above range, experience in the standard model suggests 
that gauge fields have coupling constants ad ~ 10~ 3 , 
while forces arising from scalar exchange could have a 
much wider range (down to y 2 /4-7r ~ 10~ 12 ). Thus, it is 
instructive to consider whether we have any additional 
constraints from cosmology. One possible assumption is 
that the dark matter is a thermal relic. This is a strong 
assumption, with compelling alternative models of asym- 
metric dark matter that relate the dark matter density 
to the baryon asymmetry (see, e.g. Refs. [5"TH5"B"] ). but it 
is nevertheless a general category worth examining. 

The basic requirement of a thermal relic is that 
the average annihilation cross section (<7 an n^) ~ 3 x 
10 _26 cm 3 s _1 at decoupling. Annihilation into new 
bosons is characterized by a rate constant (cr ann i>) ~ 
7ra 2 /m 2 , but it is important to recognize that there 
can be many additional annihilation channels into other 
charged states or even into other force carriers that 
simply happen not to be the dominant contributors to 
present-day scattering. Thus, this should be properly 
taken as a upper bound on ad- Thus, all but the largest 
values of ad are acceptable from the perspective of ther- 
mal freezeout, which would require a primordial asym- 
metry. 

It is important to emphasize that for our purposes, 
it is the lightest mediator that will set the character- 
istic scattering scale, not necessarily the one with the 



largest coupling constant, simply because of the satura- 
tion that sets in at low velocities. Thus, it could well be 
that the underlying Lagrangian is described by a gauge 
theory, spontaneously broken by a Higgs field, and the 
light Higgs through a small Yukawa <j) could mediate the 
self-interaction, even if the gauge interaction is mediated 
by a field too heavy to induce CDM scatterings. 

For the exothermic scenario, the parameter range for 
dark matter mass is much narrower. We can determine 
this range by requiring that the scattering in the early 
universe had ceased by the time the CDM had reached a 
low enough temperature to thermally deplete the excited 
states through the same size cross section that is relevant 
today, i.e. 

P«M< 3H , ( 7 ) 

after the time at which v x < t?dwarf ~ 10 kms -1 
(where H is the Hubble parameter). The temperature 
T x at which the WIMPs slow down to the speed in 
current dwarf galaxies is determined by the tempera- 
ture Tdoc at which they decouple from the photon bath, 
T x ~ T 2 /T dcc , where T 7 - T (l + z). Thus, set- 
ting p x — rt x p c0 (T 1 /T ) 3 , where p cQ is the present-day 

critical density and H 2 ~ (8nG/3)g*T* , yields m x ~ 
m e x MeV/Tdcc- Since it would be unnatural to have 
Idee ^ m e (as we assume the WIMPs are not charged), 
this suggests WIMPs in the sub- MeV mass range. In such 
a case, it would be most natural to have an asymmetric 
model such as described in Ref. [56 . 

Another possibility would be that the excited states are 
produced non-thermally, i.e., a heavier state freezes out 
and populates the excited states after they are adequately 
dilute that they cannot down-scatter on themselves until 
the present day. In this case, our expectation would be 
for characteristically weak-scale masses in order to have 
the appropriate relic abundance. 

From a model building perspective, such scenarios 
would most naturally be constructed within non-Abelian 
gauge theories. This is because for Abelian gauge theo- 
ries, the scattering always changes the number of excited 
states by a multiple of two, i.e., 11 •<=>■ 22 and 12 21. 
This suggests that it will be difficult to inelastically scat- 
ter the 1 states. However, in the case of a non-Abelian 
theory, more states, and thus more scattering possibili- 
ties, are present. For instance, 13 => 22 would allow an 
exothermic scattering for even ground-state particles. 

In summary, there is a significant but nonetheless con- 
strained range of parameters which could yield the YIDM 
and YIDM* phenomenology. We leave detailed model 
building for future work. 
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